Measurement of the damping of nuclear shell effect in the doubly magic 208 Pb region 
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The damping of the nuclear shell effect with excitation energy has been measured through an 
analysis of the neutron spectra following the triton transfer in the 7 Li induced reaction on 205 Tl. 
The measured neutron spectra demonstrate the expected large shell correction energy for the nuclei 
in the vicinity of doubly magic 208 Pb and a small value for 184 W. A quantitative extraction of the 
allowed values of the damping parameter 7, along with those for the asymptotic nuclear level density 
parameter a, has been made for the first time. 
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The shell effect is a cornerstone of the mean field theory 
describing finite fermionic systems. The shell structure in 
atoms decides the chemical properties of the correspond- 
ing elements. In nuclear physics the spin orbit coupling, 
in addition, plays a dominant role in deciding the shell 
closures and the associated magic numbers of protons 
and neutrons. The nuclei having such numbers of neu- 
trons and protons have an extra stability with respect 
to that expected from the average behaviour described 
by the liquid drop model (LDM). Many important nu- 
clear phenomena such as the occurrence of super heavy 
elements [TJ [2] , fission isomers [3j H] , super-deformed nu- 
clei [5] and new magic numbers in exotic nuclei [6] [7] are 
the consequences of the shell effect. The shell effect also 
affects another fundamental property of the nucleus viz. 
the nuclear level density (NLD). The NLD is an indis- 
pensable input to the statistical calculation of compound 
nuclear decay and thus an important physical quantity 
for many practical applications, such as the calculations 
of reaction rates relevant to nuclear astrophysics, nuclear 
reactors and spallation neutron sources. 

The NLD was first calculated by Bethe using a non- 
interacting Fermi gas model, without shell effects, arriv- 
ing at its leading dependence on excitation energy (Ex) 
and angular momentum (J)0[9]. The generic behaviour 
with respect to Ex is described by e 2 ^ a£;x . Here V is 
the NLD parameter which is related to the single particle 
density at the Fermi energy. Direct measurements of the 
NLD are based on the study of slow neutron resonances, 
which are mainly s- and p-wave, and are extrapolated 
to higher J values to estimate the angular momentum 
summed or total NLD. The total NLD inferred from such 
a measurement shows that on the average the level den- 
sity parameter a increases linearly with the mass number 
(A) of the nucleus as a ^A/8 MeV -1 . However, there is a 
significant departure from this liquid drop value at shell 
closures. This departure is the largest for the doubly 
magic nucleus 208 Pb, where a (at Ex ~7 MeV) is as low 
as A/26 MeV -1 . This shell effect on the NLD parameter 



is expected to damp with excitation energy so that a ap- 
proaches its liquid drop value at Ex ~ 40 MeV [10] . It is 
important to make measurements on the damping of the 
shell effect over a wide Ex range. To our knowledge, no 
such measurement has been reported. 

Experimental information on the damping of the shell 
effect can be obtained by measuring the Ex dependence 
of the NLD over a wide range, typically ~ 5-40 MeV. One 
method, which is limited to the particle bound excitation 
energy region, involves the measurement of continuum 
7-ray spectra following inelastic scattering and transfer 
reactions pTj. Both NLD and 7-ray strength function are 
inputs to the analysis of the spectra. Syed et al. used 
3 He induced inelastic scattering and single nucleon trans- 
fer reaction to populate 205 - 208 pb [12] and extracted the 
energy dependence of NLD from the coincident 7 spec- 
trum up to Ex ~6 MeV. Another method of addressing 
the Ex dependence of NLD over a wider range is by 
measuring particle evaporation spectra following heavy 
ion fusion reaction and using a statistical model analy- 
sis [13]. Lunardon et al.[I3] measured proton evaporation 
spectra in 10,11 B+ 198 Pt reactions and extracted the NLD 
in 208 Pb at an excitation energy ~ 50 MeV. However, this 
excitation energy is above the region influenced by the 
shell effect and the extracted NLD showed the expected 
liquid drop behaviour. It is indeed difficult to access a 
much lower excitation energy region using such heavy ion 
fusion reaction because of the large Coulomb barrier in 
the entrance channel. One way out of this difficulty is 
to measure particle evaporation spectra following trans- 
fer induced fusion process populating particle unbound 
states. 

In this letter, we present an exclusive measurement of 
neutron spectra from 208 Pb, following triton transfer in 
the 7 Li+ 205 Tl reaction, in coincidence with ejectile alpha 
particles. The nucleus 208 Pb (formed in the excitation 
energy range 19-23 MeV in this work) decays predom- 
inately by first step neutron emission populating 207 Pb 
in the E x ~ 3 - 14 MeV. Over this E x range, the NLD 
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FIG. 1: (Color online) Two dimensional plot of ZCT vs en- 
ergy deposited in one of the CsI(Tl) detectors in 7 Li+ 205 T1 
reaction. The inset shows the projected spectrum of the al- 
pha particles defined by the dotted two dimensional gate. The 
vertical lines define three alpha energy bins (see text). 



parameter is expected to show a significant change due 
to the damping of the shell effect. We have also made the 
above measurement with a 181 Ta target populating nu- 
clei in the 184 W region where the shell effect is expected 
to be small. 

The experiment was performed at the Mumbai Pel- 
letron Linac Facility (PLF) using a 30 MeV pulsed 7 Li 
beam of width ~1.5 ns (FWHM) and period ~ 107 ns. 
Self-supporting foils of 4.7 mg/cm 2 205 T1 (enriched to 
>99%) and 3.7 mg/cm 2 181 Ta (-100% natural abun- 
dance) were used as targets. Alpha particles were de- 
tected at backward angles (~126°-150°) in 8 CsI(Tl) de- 
tectors of dimensions 2.5 cm x 2.5 cm x 1 cm (thick) 
coupled to Si(PIN) photodiodes and placed at a distance 
of — 5 cm from the target. The detectors were covered 
with an aluminised mylar foil of thickness ~1 /im. Parti- 
cle identification was done using the standard pulse shape 
discrimination method by measuring the zero cross over 
timing (ZCT) of the amplified bipolar pulse. 

Neutrons were detected using an array of 15 plastic de- 
tectors each of dimension 6 cm x 6 cm x 100 cm viewed 
by two photomultipliers (PMT), one at each end [T5] . 
The array was placed at a mean angle of 90° to the beam 
direction and at a distance of 1 m from the target. The 
neutron energy was measured using the time of flight 
(TOF) technique. The data were collected in an event 
by event mode using a CAMAC based data acquisition 
system. The parameters recorded were (a) left and right 
timing of each plastic detector with respect to RF from 
the beam pulsing system using time to digital converters, 
(b) integrated charge of anode pulses (which is related to 
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FIG. 2: Time of flight spectrum in 7 Li + 205 T1 reaction for 
the central energy bin of alpha particles 



the electron equivalent energy, E ee , deposited in the plas- 
tic detector) from the left and right PMTs using charge 
to digital converters, (c) timing of CsI(Tl) detectors with 
respect to the pulsed beam, (d) energy deposited in the 
CsI(Tl) detectors (E Cs i) and (e) ZCT of the Csl detec- 
tors. 

A typical ZCT - Ec s i 2D-spectrum is shown in Fig. [I] 
displaying a clean separation of various groups of par- 
ticles. The energy calibration of the CsI(Tl) detec- 
tors, in the range E a ~ 5 - 25 MeV, was done using 
a 229 Th alpha source and the 12 C( 12 C, a) 20 Ne reaction 
at E( 12 C) = 24 MeV populating discrete states in 20 Ne. 
For the latter measurement, the carbon targets, backed 
by 1 - 3 mil thick Ta foils, were placed 23 cm upstream 
of the centre of the reaction chamber and the detectors 
brought to 0° to reduce the kinematic energy spread. The 
projected alpha energy spectrum for the 205 Tl target is 
shown in the inset of Fig. [I] The calibration of the en- 
ergy deposited in the plastic detector (in E ee ) was done 
using Compton tagged recoil electrons from 137 Cs and 
60 Co 7-ray sources. The time calibration was done using 
a precision time calibrator. The TOF, position informa- 
tion and geometric mean of the energy deposited for the 
neutron events in the plastic detector have been derived 
as in Ref. [15]. In order to minimize the contribution 
of scattered neutrons a TOF dependent energy thresh- 
old (increasing with decreasing TOF) was used to obtain 
the final TOF spectra after subtracting the contribution 
from the random coincidences. A typical TOF spectrum 
is shown in Fig. [2] The efficiency of the plastic detector as 
a function of incident neutron energy and energy thresh- 
old was calculated using a Monte Carlo simulation code 
[T5] . The efficiency corrected energy spectra of neutrons 
were derived from the TOF data. 

The neutron energy spectra for the Tl target are shown 
in Fig.|3ja) for three alpha energy bins, defined in Fig. [I] 
An overall decrease in the slope of the spectra with the 
increase in alpha energy (implying a decrease of Ex in 
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FIG. 3: (Color online) (a) Measured neutron spectra from 
205 T1 target for three alpha bins I, II, III (see Fig. 0, (b) 
and (c) show measured neutron spectra from 205 Tl and 181 Ta 
targets and statistical model (SM) calculations for the central 
alpha energy bin which corresponds to an excitation energy 
- 21 MeV. 



208 Pb) is consistent with the statistical nature of the neu- 
tron decay from an equilibriated nucleus. Fig. |3jb) and 
(c) show the energy spectra for both the targets gated 
by the central alpha energy bin. The excitation energy 
dependence of the NLD parameter a, which includes the 
shell effect and its damping, has been parameterised by 
Ignatyuk [16] as 
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Here a is the asymptotic value of the NLD parameter 
in the liquid drop region, A s is the shell correction en- 
ergy, which is the difference between the experimental 
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FIG. 4: (Color online) Comparison of data with SM calcu- 
lation (a) for different values of As and 7 = 0.055 MeV -1 
and (b) for different values of 7 with As = 13 MeV. The 
asymptotic level density parameter a = A/8.5 MeV -1 . 



binding energy and that calculated from the LDM , U 
is the energy available for intrinsic excitation after sub- 
tracting rotational and pairing energies from Ex and 7 
is the damping parameter. A statistical model (SM) 
calculation was done using the code CASCADE [17] 
with the above prescription for a. Figs. [3^b) and (c) 
show the calculated spectra using a = A/8.5 MeV -1 
and 7 = 0.055 MeV -1 [18 . It is seen from the fig- 
ure that a shell correction energy As = 13 MeV fits 
the shape of neutron spectrum for the Tl target while 
As = 2.2 MeV does not. An opposite behaviour is seen 
for the Ta target. These values agree with those obtained 
from the experimental nuclear masses and the calculated 
LDM values [19]. The present data, therefore, is consis- 
tent with the shell correction energies derived from the 
nuclear masses. 

We now address the extraction of damping parame- 
ter from the present data. It may be pointed out that 
constraining all three parameters, a, As and 7, is not 
possible from the data addressing even a much wider ex- 
citation energy range. By fixing any two parameters the 
third one can be constrained. For example, Fig. [4ja) 
shows the sensitivity to As for 7 = 0.055 MeV -1 and 
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FIG. 5: Exclusion plot of 5 a - 7 for A s— 13 MeV, where 
5a — A/a. The acceptable values are within the contour. 

a = A/8.5 MeV -1 . Similarly, Fig. [4^b) shows the sen- 
sitivity to 7 for A S =13 MeV and a = A/8.5 MeV -1 . 
Since the shell correction energy is known with a reason- 
ably good accuracy (within a few hundred keV [T9]). we 
have fixed As = 13 MeV and searched for an acceptable 
range of a and 7. 

The calculations were performed with 5a (= A/a) and 
7 ranging from 6.5 - 11.0 MeV and 0.02 - 0.08 MeV" 1 , 
respectively. Fig. [5] shows a 5a - 7 two dimensional ex- 
clusion plot, the region inside the contour representing 
the acceptable range of parameter values for fitting the 
present data. It can be seen from the figure that the ac- 
ceptable range of 5a lies between 8.0 and 9.5 MeV. The 
parameter 7 controlling the damping of the shell effect 
can be constrained to (0.055 ± 0.015) MeV -1 . This is 
different from the value extracted from the neutron reso- 
nance data viz. (0.079 ± 0.007) MeV -1 [20]. This could 
be due to the differences in the angular momentum states 
sampled in the two works. Moreover, the present work 
address a specific nuclear region whereas the analysis of 
Ref. [20 is global in character. 

Finally we discuss the possible sources of uncertainties 
in the present experimental method. While the major 
contribution to the a-coincident neutron spectra is ex- 
pected to arise from triton transfer-fusion reaction, there 
are other direct processes that could contribute. The 
proton pickup and 2-neutron transfer cross sections are 
small [21] and can be ignored. A Monte Carlo calcula- 
tion of the alpha-neutron coincidence spectrum reveals 
that the contribution from the one neutron and one pro- 
ton transfer is a small fraction (< 5%) in the region of 
interest, even if the cross sections are the same as that 
of the main reaction. The most relevant reaction is the 
deuteron transfer followed by 5 He breakup. However, 



the spectroscopic factor for the d+ 5 He configuration is 
expected to be much smaller than the t+ 4 He configura- 
tion [22 leading to a small contribution from this reac- 
tion. 

In conclusion, we have for the first time measured the 
effect of the shell correction on the level density param- 
eter over a range of excitation energy where the effect 
of damping is significant. The experimental results show 
that the shell correction is indeed necessary to explain 
the data and is pronounced in the Pb region. The shell 
damping factor 7 = (0.055 ± 0.015) MeV -1 has been 
extracted from the data. A precise measurement of the 
damping parameter in heavy magic nuclei will be an use- 
ful input in the current research on the formation of super 
heavy nuclei from heavy ion fusion reactions. The preci- 
sion of the present method can be improved by using a 
sharper time profile of the pulsed beam, pulse shape dis- 
crimination based neutron detectors and a larger angular 
coverage. 
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